Motional Spin Relaxation in Large Electric Fields 
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We discuss the precession of spin-polarized Ultra Cold Neutrons (UCN) and "^He atoms in uniform 
and static magnetic and electric fields and calculate the spin relaxation effects from motional i; x 
E magnetic fields. Particle motion in an electric field creates a motional v x E magnetic field, 
which when combined with collisions, produces variations of the total magnetic field and results in 
spin relaxation of neutron and ^He samples. The spin relaxation times Ti (longitudinal) and T2 
(transverse) of spin-polarized UCN and '^He atoms are important considerations in a new search 
for the neutron Electric Dipole Moment at the SNS nEDM experiment. We use a Monte Carlo 
approach to simulate the relaxation of spins due to the motional v x E field for UCN and for ^He 
atoms at temperatures below 600 mK. We find the relaxation times for the neutron due to the v x E 
effect to be long compared to the neutron lifetime, while the ^He relaxation times may be important 
for the nEDM experiment. 



I. INTRODUCTION 

The phenomenon of longitudinal and transverse spin re- 
laxation due to magnetic field inhomogeneities has been 
discussed extensively [H 0, S, 0, @|- There, the fiuctuat- 
ing magnetic fields experienced by the particle spins dur- 
ing transport through magnetic field gradients coupled 
to collisions with other particles or a holding cell's walls 
lead to relaxation. For a combined magnetic and electric 
field configuration, another possible source of spin relax- 
ation has been considered by Lamoreaux [1| and results 
from fluctuations in the motional magnetic fleld in the 
particle's rest frame. Here we present a more detailed 
discussion of this effect. 

Transformation into the rest frame of a particle moving 
in an electric fleld gives rise to a motional magnetic fleld 



B^ = Ex vjc? 



(1) 



The motional magnetic fleld is always perpendicular to 
both the direction of motion and the electric fleld. For 
parallel electric and magnetic field configurations (as are 
typically employed in particle electric dipole moment 
searches) , the superposition of and a magnetic hold- 
ing field So leads to a small tilt in the total magnetic 
field B = Ba -I- B^ (see Fig. [TJ, as the motional mag- 
netic fields are typically very small relative to i?o- For 
example, in an upcoming search (n£^£'M)[3| for the neu- 
tron electric dipole moment at the Spallation Neutron 
Source (SNS) in the Oak Ridge National Laboratory, 
the motional magnetic field experienced by the neutrons 
from the 50kV/cm electric field will be on the order of 
5.6 X 10~'*mGauss per m/s of velocity, whereas the mag- 
netic holding field Bq will be about 10 mGauss. 

Nevertheless, collisions result in an abrupt change in 
the velocity vector, causing the motional magnetic field 
B„ to change both in direction and magnitude. Although 
the changes in orientation of the total magnetic field fol- 
lowing a single collision are rather small (as i?o), 
after each collision, a particle's spin will precess about 
a sHghtly reoriented total magnetic field. After a large 
number of such collisions, each of these small "kicks" in 



the field will contribute to the overall global relaxation 
of a spin ensemble. 



In the remainder of this paper, we calculate the mag- 
nitude of spin relaxation due to the v x E effect, es- 
pecially as relevant for the conditions of the upcoming 
nEDM search, which will utiHze Ultra Cold Neutrons 
(UCN, neutrons with speeds < 7 m/s) and an in-situ po- 
larized ^He "co-magnetometer" immersed in a superfluid 
^He bath at temperatures in the range 100 — 600 mK. 
Spin-relaxation times signiflcantly longer than the neu- 
tron lifetime (limiting the measurement time) are desired. 



The results in this paper are complementary to 
the studies on spin relaxation due to magnetic fleld 
gradients[ll0, S;[3| and calculations of the linear elec- 
tric fleld geometric phase effects [1, [3, both of which 
can be important considerations for the measurement of 
EDMs. 
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Figure 1: (a) The total magnetic field B as the superposition 
of the holding field Bo and the v x E motional field . (b) 
Rotation of motional field B^ due to wall collision. 
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II. THEORY 



He Atoms 



We have followed the formalism of McGregor (sl to 
describe the relaxation process of the v x E motional 
field. This uses the field auto-correlation function to cal- 
culate longitudinal (Ti) and transverse (T2) spin relax- 
ation times. The spin relaxation times for particle spins 
due to a changing B field can be calculated using the 
field auto-correlation function. With this method[5||, the 
longitudinal spin relaxation rate can be expressed as 



1 



7 



[Sbx{(^o) + SByiuJo)] 



(2) 



where 7 is the gyromagnetic ratio and Sbx{^^o) and 
Ssyii-^o) are the Fourier transforms of the field auto- 
correlation function. For a perturbing field B, the fre- 
quency spectrum is given by 



Sb, (w) 



B,[t)B,{t 



(3) 



for each of the fiuctuating field components. The average 
inside the integral in Eq. [3] is calculated for particles in 
the measurement cell colliding with the walls of the cell 
and with excitations in the superfiuid ''He. 

For the case of particles in a uniform magnetic hold- 
ing field in the z direction and moving at velocity v in 
a uniform electric field E parallel to Bo, the non-zero 
perturbation magnetic field terms are 



Bj, -Evyjf? 
By = Ev^/c^ 



(4) 



The field auto-correlation function, for this case, is pro- 
portional to the velocity auto-correlation function, as 



'B,{t)B.,,{t + T)\ = {vy{t)vy{t + T)) Eye" 
By{t) By{t + t)) = {v,{t)v,{t + T)) Eye" 



(5) 



For this work we adopt the velocity auto-correlation 
functional 



(6) 



where Tc is the mean time between collisions. Eq. [6] is 
applicable in determining the frequency spectrum when 
the Larmor frequency loq is larger than 1/td, where td 
is the diffusion time. This is the case for ^He atoms at 
temperatures T > 400 mK, where, as discussed in Sec. 
IIII A| coUisions with excitations in the superfiuid are fre- 
quent and atoms take a relatively long time to travel to 
the walls of the holding cell. The corresponding field 
auto-correlation frequency spectrum is 



Sbx{uj) 



1 + cjV2 c" 



(7) 



and a similar formula for SByi^^)- With Eq. [21 we find 
the longitudinal spin relaxation time to be 



3 1 



2 72 E^v^^^.T, 



(1 + ^0 ) 



(8) 



where we have used (w^^ + (f 2) = 2 vfj^^/S. As discussed 
in Sec. IIII A^ the mean time between collisions drops dra- 
matically as the temperature of the '^He sample increases, 
and, in the nEDM setting, WqTc <C 1 for T > 300 mK, so 
that Eq. [8] can be approximated to be 



Ti 



3 1 



2 72 £;2 y2 _ 



(9) 



leaving the relaxation times independent of the magnetic 
holding field Bq for this regime. 

The transverse spin relaxation rate can be expressed 
asfl 



1 



1 

2J\ 



(10) 



but, for the case ofv x E fields in which the perturbation 
field has no z-component, the term Sbz{0) = 0, so we 
find that T2 = 2Ti. 



B. Ultra Cold Neutrons 

Although Eq. [2] can be applied to any perturbative 
field, as long as the proper expression for the auto- 
correlation function is used, we note that the earlier for- 
maHsm detailed in GambHn and Carver [l| is useful for 
the special case of UCN whose velocities do not change 
in the nEDM measurement cells. Following their proce- 
dure to calculate the spin relaxation in inhomogeneous 
fields, which does not rely on the field auto-correlation 
function, the expression for the longitudinal relaxation 
time Ti due to motional v x E field can be calculated to 
be 



3 (1.19)2AS2c'* l-H(woTc)^ 



^3 ^2 



(11) 



This formula is equivalent to Eq. [HI apart from small 
constant terms. It is also interesting to note that, for the 
case of UCN in the nEDM experiment, the mean free 
path A in Eq. [iT] is a constant of the geometry of the 
measurement cell, since UCN do not interact with other 
particles. In addition, wqTc > 1 for UCN speeds, so that 
the last factor in Eq. [lT]is roughly constant and of order 
unity. In this case, Eq. [lT]can be approximated to be 



3 (1.19)^ XBjc" 

4 ' „3 £;2 



(12) 



emphasizing the Ti oc relation. The relaxation time 
can be lengthened by increasing the magnitude of the 
magnetic holding field Bq. However, as seen in Sec. IIIIB| 
the motional spin relaxation times for UCN in the nEDM 
experiment are long compared to the neutron lifetime. 
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Figure 2: (Color online) Simulation results of the longitudinal 
polarization Pz for '^He sample in a thermal bath at T = 
100 mK vs. time. An exponential fit to the simulation data 
is also shown. The fit yields the relaxation time Ti. 



III. SIMULATION 

We have simulated the spin relaxation of UCN and 
^He atoms due to the v x E effect for uniform and 
static magnetic and electric fields. Our studies are 
based on the precession of classical spins in a mag- 
netic field. Neutrons (and '^He atoms) are confined in- 
side the rectangular nEDM measurement cell of dimen- 
sions {Ax — 10.2cm, Ay — 50 cm, Az — 7.6cm), where 
both magnetic and electric fields are oriented in the z 
direction. The mean distance between wall collisions can 
be approximated by ~ AV/A for a rectangular mea- 
surement cell, where V and A are the volume and the 
surface area of the cell respectively. For the nEDM cell, 
\w — 8 cm, which is comparable to the mean distance 
between wall coUisions of Xsim = 7.5 cm obtained from 
the simulations. As discussed in Sec. IIII A^ '^He atoms 
can interact with excitations in the superfiuid, so that 
the mean time between coUisions Tc is determined both 
by diffusion and wall collisions. 

As we neglect the effect of gravity, changes in velocity 
are only due to collisions, hence, between collisions, the 
total magnetic field experienced in the rest frame of a 
particle stays constant. We have used an analytic solu- 
tion to the equation of motion of the spin in a constant 
magnetic field 

<T = 7(T X B (13) 

to calculate the spin precession between collisions. The 
uniform magnetic field solution to Eq. [13] can be found 



to be: 

(T{t) = [ctq - i^2J|)^) cos(wO + ^SXB sin(Ljt) + 
, (g-oB)B 

(14) 

We have studied the case of totally diffuse wall col- 
lisions for both UCN and ^He atoms. The diffuse wall 
collisions are in part responsible for the spin relaxation. 
Our simulations have shown that, for the v x E studies 
described in this paper, purely specular wall collisions, in 
the absence of particle diffusion, do not result in signifi- 
cant longitudinal relaxation of the spin polarization. In 
this case, the velocity of a particle after a wall collision 
is highly correlated to its initial velocity and the times 
between wall collisions are repetitive, so that the mo- 
tional field By changes in a repeatable way and does not, 
generally, allow the spins to precess far from their initial 
states. We were able to determine from the simulation 
data that only a small fraction (< 10~^) of particle spins 
were precessing away from their initial states for the case 
of specular wall collisions. We believe that this effect 
occurs because these particles are traveling in resonant 
paths, where the correlation of velocities and precession 
times between wall collisions is such as to maximally ro- 
tate the spins out of alignment. Furthermore, the walls 
of the measurement cells in the nEDM experiment are 
not likely to be purely specular. 

In the simulations we track the spin state of each parti- 
cle and calculate the longitudinal, Ti, and transverse, T2, 
relaxation times due to the v x E effect for single speed 
UCN and for ^He atoms at temperatures below 600 mK. 

In the longitudinal spin relaxation simulations, parti- 
cles start with their spins aligned with the holding field 
Bq. Due to the vxE effect, the ensemble of spins relaxes 
exponentially with a time constant Ti, as shown in Fig. 
[21 In contrast, in the transverse spin relaxation simula- 
tions, particles start with their spins perpendicular to the 
holding field Bo and the spins relax with an exponential 
time constant T2. 

A. ^He atoms 

We have simulated the spin relaxation due to the vxE 
effect in a magnetic holding field Bq — \Q mGauss and an 
electric field Eq — 50kV/cm, for a sample of '^He atoms 
at temperatures below 600 mK. UnHke UCN (see Sec. 
nil Bp . whose speeds were constant in our simulations, 
the '^He sample obeys a Maxwell-Boltzmann velocity dis- 
tribution of the form 

where T is the temperature of the sample, k is the 
Boltzmann constant, and ~ 2.2 M3 is the effective 
mass of ■^He atoms in superfiuid ^He, with M3 as the 
■^He atomic mass[llj. The mean speed of '^He atoms is 
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Figure 3: (Color online) The spin relaxation times Ti and T2 
for ■'He atoms as a function of temperature, plotted together 
with the theoretical formulae derived from Eq. O 
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Table I: Temperature dependence of velocity and mean free 
path for a sample of '^He atoms at selected temperatures. The 
4mK point represents ^He velocities comparable to UCN. 



V = (v) — y/8 kT/n A'h, while the root mean square 
speed is Vrms = ^/W) = v^3fcT/Af*. 

The colUsions with the walls and interactions with the 
superfluid were implemented to be totally diffuse. In ad- 
dition, we included thermalization of '^He atoms, so that 
the magnitude of the velocity is randomly sampled from 
the original velocity distribution (Eq. [15]) after each colli- 
sion with either phonons or the walls of the measurement 
cell. 

The '^He atoms can interact with excitations in the su- 
perfluid "^He in the nEDM cell. Due to the low fractional 
density of ^He atoms, typically on the order of one '^He 
atom per 10^" atoms of *He, collisions with other '^He 
atoms are negligible. The mass diffusion coefficient for 
■^He atoms in superfluid He depends sensitively on the 
sample temperature and was found experimentallvfl^ to 
be 



D 



(1.6 ±0.2) cm2 -s- 



(16) 



Since the phonon velocity in superfluid ^He is much 
higher than the velocity of individual atoms the 
mean time between phonon interactions with ^He atoms 
in the measurement cell is constant for a given temper- 
ature and independent of the instantaneous '^He velocity 
and can be expressed as[l3| 



3D 



(17) 



The effective mean time between collisions Tc, which 
takes into account both superfluid excitations and wall 
collisions, is 



1 

Tr 



1 



(18) 



The mean time td that it takes for '^He atoms to diffuse 
to the walls of the measurement cell is useful for defining 
the regime in which the formalism described by Eq. [6] is 
valid, and can be expressed as 



TD 



(19) 



where i? ~ 8 cm is the characteristic size of the measure- 
ment cell, or the radius for a spherical cell. 

Table H] shows the mean distance between interactions 
with excitations in the superfluid Ap = TpV, for selected 
simulation temperatures, along with other relevant quan- 
tities. In the low temperature regime (T < 200 mK) the 
mean distance between collisions with excitations Xp is 
longer than the dimensions of the nEDM measurement 
cell, so wall collisions dominate and the mean free path 
is the mean distance between wall collisions. In the high 
temperature regime (T > 300 mK) Ap < 5 cm, so colli- 
sions with excitations dominate. When the collision fre- 
quency is high relative to the precession frequency, the 
spin relaxation times due to the v y. E effect increase. In 
the mid-temperature regime (T ^ 100 — 300 mK), where 
the quantity ujqTe) is of order unity, Eq. [7] is not a good 
approximation for the frequency spectrum, and this is 
where we see the largest deviation between the simula- 
tions and the predicted relaxation times. For u)qTd ^ 1, 
a full diffusion theory calculation is necessary to deter- 
mine the frequency spectrum. In any case, the nEDM 
geometry and fields configuration is such that the quan- 
tity ijJqTu ^ 1 for all ranges of temperatures in our study, 
keeping the deviation small. 

We analyzed the data from our simulations to calcu- 
late the spin relaxation times Ti and T2 for ^He atoms 
for temperatures below 600 mK. The results are shown in 
Fig. [21 where relaxation times are plotted as a function of 
temperature. The simulation data agree with the theo- 
retical diffusion calculations, especially for temperatures 
higher than 400 mK, in which the collision frequency 1/tc 
is large and the diffusion times tjj are long. The simula- 
tion data also agree with the relation T2 = 2Ti. 

Our calculations and simulations have shown that the 
longitudinal spin relaxation times Ti for '^He atoms at 
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Figure 4: (Color online) Simulation results for longitudinal 
spin relaxation times Ti for UCN as a function of velocity 
and with different electric field strengths. The simulation 
data agree with the theoretical formalism described by Eq. 

EU 

temperatures in the range 100 — 400 mK can be as short 
as 10'' s. These relaxation times are comparable to the 
measurement time of nEDM and may be an important 
component of the total relaxation rate. In this regard, 
the temperature is a parameter of the measurement that 
can be chosen to optimize spin relaxation times. 

B. Ultra Cold Neutrons 

We have simulated the spin relaxation due to the vxE 
effect for UCN for the conditions in the nEDM experi- 
ment. The measurement cell, as designed for the nEDM 
experiment, has a UCN Fermi potential of 134 neV, which 
is equivalent to a maximum UCN velocity of 3.6 m/s and 
constitutes an upper limit for the velocity spectrum of 
UCN. 

The relaxation times are found to be shorter for higher 
UCN speeds, as discussed in Sec. Ill Bl We have sim- 
ulated UCN at different velocities (see Fig. IH, ex- 



periencing diffuse collisions with the walls of the cell. 
The UCN do not interact with other particles in the 
cell so their speed stays constant in our simulations. 
The nEDM magnetic field configuration was used, with 
Bq — lOmGauss. 

The results from the UCN simulations agree with the 
theoretical results (Eq. [TT|) in the regime in which colli- 
sion frequencies are low with respect to the Larmor fre- 
quency. The expected dependence of Ti for UCN 
can be seen in Fig. [H for two values of -B = 5kV ■ cm~^ 
and 50 kV • cm~^. The relevant relaxation times for UCN 
with speeds less than 4 m/s are on the order of 10^ s, for 
the configuration of the nEDM measurement. These re- 
laxation times are much longer than the neutron lifetime 
and should not constitute a limiting factor for the exper- 
iment. 



IV. CONCLUSION 

We have studied the longitudinal and transverse spin- 
relaxation times Ti and T2 from motional v x E fields. 
The results from Monte Carlo simulations and theo- 
retical calculations of the relaxation times are shown 
to be in agreement. We investigated this phenomenon 
within the context of the operating parameters for an 
upcoming search for the neutron electric dipole moment 
utilizing Ultra Cold Neutrons and a polarized '^He co- 
magnetometer. For UCN, the relaxation times were 
found to be on the order of 10^ s, much longer than 
the neutron lifetime. In contrast, our calculations have 
shown that the relaxation times for "^He atoms will vary 
strongly with temperature, and can be as low as lO'^s. 
Thus, consideration of the relaxation effects from v x E 
motional magnetic fields can infiuence the choice of the 
experiment's operating temperature. 
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